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The (n, 7) cross section of ^"^Pb has been measured at the CERN n_TOF facihty with high 
resolution in the energy range from 1 eV to 620 keV by using two optimized CeDg detectors. 
In the investigated energy interval about 130 resonances could be observed, from which 61 had 
enough statistics to be reliably analyzed via the R-matrix analysis code SAMMY. Experimental 
uncertainties were minimized, in particular with respect to (i) angular distribution effects of the 
prompt capture 7-rays, and to (ii) the TOF-dependent background due to sample-scattered neutrons. 
Other background components were addressed by background measurements with an enriched '^'^^Pb 
sample. The effect of the lower energy cutoff in the pulse height spectra of the CeDg detectors was 
carefully corrected via Monte Carlo simulations. Compared to previous ^"^Pb values, the Maxwellian 
averaged capture cross sections derived from these data are about 20% and 9% lower at thermal 
energies of 5 keV and 30 keV, respectively. These new results have a direct impact on the s-process 
abundance of ^"^Pb, which represents an important test for the interpretation of the cosmic clock 
based on the decay of ^^*U. 

PACS numbers: 25.40.Lw,27.80.+w,97.10.Cv 

Keywords: Neutron capture cross sections; Nuclear astrophysics; Pulse height weighting technique; CeDg 
scintillation detectors; Monte Carlo simulations 



I. INTRODUCTION 

Similar to the majority of the stable isotopes beyond 
iron, 206,207,208p|-, g^^^ 209gj g^j.g synthesized by the rapid 

(r-) and slow (s-) neutron capture processes. However, 
this mass region is particularly interesting because the 
r-process abundances are dominated by the decay of the 
short lived a-unstable transbismuth isotopes pj. This 
feature provides an important consistency check for the 
r-process abundance calculations in the actinide region, 
since the integrated r residuals are constrained by the 
difference between the solar abundance values and the 
respective s-process components. Reliable r-process cal- 
culations are required for the interpretation of the ob- 
served Th and U abundances in the ultra metal-poor 
(UMP) stars of the Galactic halo. Since these stars are 
considered to be as old as the Galaxy, the observed Th 
and U abundances can be used as cosmo-chronometers, 
provided the original Th and U abundances are inferred 
from r-process models. This dating mechanism has the 
advantage of being independent of the yet uncertain r- 
process site P, [H, HI . 

Apart from its relevance for establishing the basic 
constraints for the r-process chronometry in general, 
^"^Pb contains also dating information in itself. The 
206p|^^238u cosmochronometer was first introduced by 
Clayton in 1964 The ^^^U produced by the r pro- 
cess decays with a half life of ti/2 = 4.5 x 10^ yr over a 
chain of a and (3 decays ending at ^'^^Pb. Therefore, its 
radiogenic abundance component, N"^^^ , can be used to 
constrain the age of the parent isotope ^■^^U, and hence 
the age (A^) of the r-process. Unlike the more direct r- 
process abundance predictions derived from the Th and 
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U abundances in UMP stars, this procedure requires a 
Galactic evolution model, which describes the supernova 
rate or the frequency of the r-process events [5|. The 
drawback of this clock arises from the difficulty to iso- 
late the cosmo-radiogenic component of ^°^Pb accurately 
enough from the additional abundance components. 

Apart from these astrophysical aspects, the neutron 
capture cross section of ^°^Pb is also of importance for the 
design of fast reactor systems based on a Pb/Bi spallation 
source. Because 24.1% of natural lead consists of ^°^Pb, 
its {n, 7) cross section influences the neutron balance of 
the reactor @. 

There have been several measurements of the 
^''^Pb(n, 7) cross section, which show discrepancies that 
are difficult to understand (see Sec. IIV Ap . The aim of 
this work is to perform a new independent measurement 
with higher accuracy and in this way to determine the s- 
process contribution to the ^''^Pb abundance, A^s"®, more 
reliably. 

In fact, the s-process abundance of this isotope is al- 
most completely determined by the stellar (n, 7) cross 
section, nearly independent of the stellar model used 0] ■ 
Therefore, the uncertainty of Nf^^ arises mostly from the 
cross section uncertainty. 

Potential sources of systematic error have been sub- 
stantially reduced in the present measurement, which was 
performed at the CERN n_TOF installation. The new 
setup, and in particular the detectors themselves, were 
optimized for very low neutron sensitivity. Furthermore, 
the detectors were mounted at ^^125° with respect to the 
incident neutron beam in order to minimize the correc- 
tion for angular distribution effects. The experimental 
details are presented in Sec. HH followed by the adopted 
data analysis procedures and an evaluation of the var- 
ious systematic uncertainties in Sec. IIIII The deduced 
resonance parameters and the corresponding Maxwellian 
averaged capture cross sections in the stellar temperature 
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regime are presented in Sec. lIVl Based on these new data, 
first astrophysical implications for the s-process abun- 
dance of ^"^Pb are discussed in Sec. IVl 



II. MEASUREMENT 

The time-of-flight (TOF) measurement was performed 
at the CERN n_TOF instaUation Q using a set of two 
CeDg detectors. Neutrons were produced by a 20 GeV 
proton beam on a lead spallation target. The spalla- 
tion source was surrounded by a 6 cm thick water layer, 
which served as a coolant and as a moderator for the 
initially fast neutron spectrum. The beam was char- 
acterized by intense bunches of (3 to 7)xl0^^ protons, 
a width of 6 ns (rms), and a repetition rate of only 
0.4 Hz. This extremely low duty-cycle allows one to 
perform (n, 7) measurements over a broad neutron en- 
ergy interval from 1 eV up to 1 MeV and to achieve 
favorable background conditions. Data were recorded by 
means of an advanced acquisition system with zero dead 
time, based on 8 bit Flash-Analog-to-Digital Converters 
(FADC), with 500 MHz sampling rate and 8 MB buffer 
memory Q. 

The measurement was performed with an enriched 
metal sample 8.123 g in mass and 20 mm in diameter. 
The sample was enriched to 99.76% in ^"^Pb with small 
impurities of ^o^pb (0.21%) and ^ospb (0.03%). 

Capture events were registered with two CgDg 7 -ra.y 
detectors optimized for very low neutron sensitivity [lOj. 
A sketch of the experimental setup is shown in Fig. 2 of 
Ref. [ill] . The absolute value of the neutron fluence was 
determined by regular calibration measurements with an 
0.5 mm thick gold s amp le and by using the saturated 
resonance technique [12] for the first gold resonance at 
En — 4.9 eV. The energy differential neutron flux was de- 
termined with a relative uncertainty of ±2% from the flux 
measurement with a 235,238jj gggjon chamber calibrated 
by Physikalisch-Technische Bundesanstalt (PTB) p^ . 
The neutron intensity at the sample position was also 
monitored by means of a 200-/ig/cm^ thick ^Li foil in the 
neutron beam about 2.5 m upstream of the capture sam- 
ple. The ^Li foil was surrounded by four silicon detectors 
outside of the beam for recording the '^H and a particles 
from the (n, a) reactions. 

Compared to previous measurements [T3 . [Tsj . the 
present setup had the advantage that the detectors were 
placed at ~125° with respect to the incident neutron 
beam. In this way, the corrections for angular distribu- 
tion effects of the prompt capture 7-rays were strongly 
reduced. This configuration led also to a substantial re- 
duction of the background from in-beam 7-rays scattered 
in the sample [l^. 



III. CAPTURE DATA ANALYSIS 

The response function of the CeDg detectors needs to 
be modified such that the detection probability for cap- 
ture cascades becomes independent of the cascade mul- 
tipolarity. This was accomplished by application of the 
Pulse Height Weighting Technique (PHWT) Based 
on previous experience [Til . [Tsl . Il9l | , the weighting func- 
tions (WFs) for the gold and lead samples were obtained 
by means of Monte Carlo calculations. The accuracy 
of the WFs was verified with the method described in 
Ref. [1^, by which the calculated WFs were applied to 
Monte Carlo simulated capture 7-ray spectra. Using this 
procedure, the uncertainty of the WFs was estimated to 
be smaller than 0.5% for the samples used in the present 
experiment. 

The weighted count rate iV™ is then transformed into 
an experimental yield. 



where the yield-normalization factor is determined 
by calibration measurements using the saturated 4.9 eV 
resonance in gold. Nn denotes the neutron flux and Ec 
the effective binding energy. 

The yield in Eq. ^ is still subject to several correc- 
tions. The common effects of the background and of the 
low energy cutoff in the pulse height spectra of the 7 de- 
tectors are described in Secs. lIII Xl and llTl Bl respectively. 
The measurement on ^°^Pb is particularly sensitive to the 
angular distribution of the prompt capture 7-rays. The 
impact of this effect is described in Sec. IIII CI 



A. Backgrounds 

A major source of background is due to in-beam 7- 
rays, predominantly from neutron captures in the water 
moderator, which travel along the neutron flight tube 
and are scattered in the ^°^Pb sample. This background 
exhibits a smooth dependence on neutron energy, with a 
broad maximum around En ~ 10 keV. The shape of this 
background was determined from the spectrum measured 
with an isotopically pure ^osp^ sample, which contains 
practically no resonances in the investigated neutron en- 
ergy range. This spectrum was properly scaled and used 
as a point-wise numerical function in the R-matrix anal- 
ysis of the ^°^Pb capture yield (see Sec. lIVp. 

Another type of background arises in the analysis of 
resonances with a dominant neutron scattering channel, 
r„ >> Tj. In such cases, there are about Tn/T^ scat- 
tered neutrons per capture event. These scattered neu- 
trons can be captured in the detectors or in surrounding 
materials, thus mimicking true capture events. This ef- 
fect was estimated to be negligible for all the resonances 
listed in Table [U 
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B. Digital threshold 

As mentioned in Sec. |TT1 FADCs were used for record- 
ing directly the analog output signals of the CgDg detec- 
tors. Without any further discrimination, 8 MB of data 
would have been acquired per proton pulse in each de- 
tector. Depending on the sample, this enormous amount 
of data could be reduced by factors of 20 to 100 by us- 
ing a zero suppression algorithm (see Ref. [9] for details). 
By this method events below a certain pulse-height are 
discriminated by a constant digital threshold analogous 
to conventional data acquisition systems, where an elec- 
tronic threshold is used to reduce backgrounds and dead 
time effects. 

Due to this threshold, the pulse height spectra of the 
CgDg detectors exhibit a low energy cutoff at a certain 
value of the signal amplitude (see Fig. [1]). In this experi- 
ment the threshold was set at a 7-ray energy of 320 keV. 
If the pulse height spectra of the ^"^Pb sample and of 
the gold sample used for normalization would have the 
same shape, the fraction of weighted counts below this 
threshold would nearly cancel out in the expression for 
the yield, 



Y^320 keV Ti/Pb pPb , V^Bc T/f/Pb pPb 

yexp ^ l^OkeV i Z^320 keV i 

E320 he.V Ti/Au dAu i V^^^c T/f/Au nAu 

OfceV I "I" Z^320 keV i 

X^Ea M/Pb pPb 

^ Z^320 keV i 

Z^320 keV ^^i 

Here, the Wi and Ri are the corresponding weighting fac- 
tors and response functions for a certain time of flight 
channel, respectively. However, this approximation is 
only valid within 4 to 5%, because the pulse height spec- 
tra of captures on 206p|-, ^^^^ ^^'^Au differ significantly 
near threshold (Fig. [J). 

This effect has been taken into account in the deter- 
mination of the experimental capture yield by simulating 
the capture cascades of each isotope as described in detail 
in Refs. [ll|, dj [S| ■ Fig. [1] shows that the experimental 
spectra above the digital threshold are well reproduced 
by the simulations. With this correction the experimen- 
tal yield becomes 




Edep (MeV) 

FIG. 1: Pulse height spectra for the 4.9 eV resonance in gold 
(grey) and for the 3.3 keV resonance in ^"^Pb (black), arbi- 
trarily scaled. The dashed lines are the MC-calculated 7-ray 
spectra for the two resonances. The linear scale used in the 
inset illustrates the large difference between the simulated 
spectra below a threshold of 300 keV. 

C. Angular distribution effects 

Neutron capture with orbital angular momentum I > 
leads to an aligned state in the compound nucleus, per- 
pendicular to the direction of the incident neutron. Given 
the small multiplicity (m = 1 to 2) of the capture cas- 
cades in ^"^^Pb, most of the prompt 7-rays registered with 
the CgDe detectors still carry this anisotropy, which af- 
fects the measured yield. The angular distribution is in 
general given by, 

W{e) = AkPkicosO) = 1 + A2P2{cos9) + 

k 

+A4,P4{cos9) + A(iPe{cos9), (4) 

where Pk{cos9) are the Legendre polynomials of order k 
and Ak are coefficients, which depend on the initial (J) 
and final (J') spin values, on the multipolarity (Z) of the 
transition, and on the degree of alignment. The angular 
distribution effects in the capture yield are minimized 
(although not avoided) by setting the detectors at 125°. 
Since each CgDg detector covers a substantial solid angle, 
capture 7-rays are registered around 125° ± A9. For the 
actual setup of the present measurement one finds A9 w 
28°. 



ft \^Ea T/f/Pb pPb 

yexp ^ ^Pb Z^320 keV i fo\ 
ft Y^Ea wAu nAu ' ^ ' 

■'Au 2^320 keV 

For the adopted digital threshold the yield of the 
4.9 eV resonance in ^^^Au needs to be scaled by a fac- 
tor /^^ — 1.071(3), whereas the yield of the resonances in 
206pj^ required a correction of /pj^ — 1.021(5) due to their 
harder spectrum. Hence, the correction factor of the final 
experimental yield was /* = /pb/Zlu = 0.952(4). 



1. Resonances with spin J = 1/2 

For resonances with J = 1/2 it can be assumed that 
they decay directly to the ground state (J'^ = 1/2") 
or to the first or second excited states with = 5/2" 
and = 3/2~, respectively (see also Fig. [2]). In these 
cases, one finds that A2 = Ai = Aq = Q. Therefore, only 
resonances with spin J > 1/2 may be affected by angular 
distribution effects. 
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TABLE I: Measured decay patterns from resonances with spin 
J = 3/2 The systematic uncertainty in the yield of 

each resonance due to the angular distribution of the involved 
transitions are given in the last column. 



E„ (keV) 




Intensity (%) 




3/2" 








(keV) 








6737.9 


6168.6 


5840.8 


4114.5 




3.36 


76.0(27) 


2.5(8) 


8.58(11) 


13.0(8) 


±10% 


3.36" 


60 


2.5 


24.5 


13 


±8% 


10.86 




100 






±2% 


21.87 




100 






±2% 


42.07 






100 




±10% 



This work. 



2. Resonances with spin J = 3/2 

In order to quantify the uncertainty due to the angu- 
lar distribution of the prompt 7-rays emitted from ex- 
cited states with J'^ = 3/2~ the de-excitation patterns 
reported in Ref. have been used (Table |T| . 



6738.2 



(d) 



5/2+ 



(c) 



3/2- 
5/2- 

1/2- 



(b) 



(a) 



2623.1 



897.5 
' 569.6 

0.0 



FIG. 2: Level scheme and decay patterns for 
energies are in keV. 



^Pb hj. AU 



of the angular distribution effects arises from the un- 
known admixtures of different multipolarities (M1+E2) 
for the transitions connecting the original excited state 
J'^ = 3/2" with any of the three lowest states (paths (a), 
(b) and (c) in Fig. [2]). As shown in Table HI the decay 
pattern and the corresponding effect on the capture yield 

CTg vary abruptly from one resonance to another. It 
is therefore difficult to assess a common systematic un- 
certainty for the remaining 3/2~ resonances. Assuming 
that the four resonances listed in Table [J constitute a 
representative sample, one may consider their standard 
deviation of cr = 4% as a realistic estimate of the system- 
atic uncertainty due to angular distribution effects. 

Resonances with J'^ — 3/2+ can be assumed to decay 
directly to the ground state through an El transition. In 
this case we have estimated an effect of 10% in the cap- 
ture yield with respect to the isotropic case. However, 
since 3/2+ resonances appear at a relatively high neu- 
tron energy, the final effect in the MACS is practically 
neghgible (see below Sec. IIVB|) . 



3. Resonances with spin J = 5/2 

For resonances in ^"^Pb with J'^ = 5/2+ the most 
probable decay would be through an electric dipole tran- 
sition to the first excited state with J'^ = 5/2~ and/or 
to the second excited state with J'^ = 3/2" (paths (b) 
and (c) in Fig. Under these assumptions, the ef- 
fect on the capture yield would be -12% for path (b) 
and 9% for path (c). However, mixtures of both decay 
paths would partly compensate the correction for angu- 
lar distribution effects. Adopting one standard devia- 

tion of the two extreme cases ~ 10% would, there- 

fore, represent a rather conservative estimate of the cor- 
responding uncertainty. Nevertheless, even such a rela- 
tively large uncertainty for the cross section of J'^ = 5/2+ 
resonances would have negligible consequences for the 
Maxwellian averaged cross section because these reso- 
nances contribute very little to the total capture cross 
section (see Sec. IIVB"]) . 



For the first resonance at 3.36 keV, fair agreement has 
been found between the relative intensities of Ref. [l3| 
and the rather coarse values deduced from the experi- 
mental pulse height spectrum (Table H] and Fig.[T]), which 
suffer from uncertainties due to background subtraction, 
limited counting statistics and poor energy resolution of 
the CeDg detectors. Therefore, an uncertainty of about 
20% has to be ascribed to the quoted 7-ray intensities. 

The estimated effect of the angular distribution on 

the capture yield {0^^ ) is given in the last column of 
Table ID These values were obtained via Monte Carlo 
simulations of the experimental setup, using the ener- 
gies and intensities listed in Table U and the prescription 
of Ref. [l^]. The main uncertainty in the calculation 



D. Summary of uncertainties 

With the WFs calculated via the Monte Carlo tech- 
nique, the accuracy of the PHWT has been investigated 
in detail by the n_TOF collaboration [18]. It has been 
shown that the capture yield can be determined from 
the measured raw data with an accuracy better than 2%. 

Other sources of systematic uncertainty pertaining to 
this measurement are due to the energy dependence of 
the neutron flux (±2%) and to the background due to 
in-beam 7-rays (±1%). In the particular case of the 
(n, 7) cross section of ^°^Pb, the uncertainty introduced 
by the angular distribution of the capture 7-rays has to 
be considered as well. This effect has been estimated to 
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contribute an uncertainty of ±4% for resonances with 
J'^ — 3/2^ and less than ±10% for resonances with 
,r = 3/2+, 5/2+. 



yield was fitted to the corrected experimental yield by 
variation of the capture width and/or neutron width 

r,,,, 



IV. RESULTS 

A total of 61 capture levels were analyzed in the neu- 
tron energy range from 3 keV up to 570 keV using the 
R-matrix code SAMMY [2l|. In the analysis, the or- 
bital angular momenta I and the resonance spins J were 
adopted from Ref. [13]. Some of the I and J parameters 
listed in Table |TT] are tentative or arbitrary if missing in 
Ref. [l^l. We list all the parameters used in our analy- 
sis so that the final values can be recalculated if neces- 
sary. The capture yield ^(iJo, r„, F^) was parameterized 
with the Reich-Moore formalism, and a channel radius of 
9.5 fm was used for all partial waves. This parameterized 



/* xr-^^f = B-f y(£;o,F„,F^), 



(5) 



where /* is the global yield correction factor given in 
Sec. IIII Bl The term B describing the background was 
parameterized as an analytical function of the neutron 
energy in the range between 1 eV and 30 keV. Beyond 
30 keV, B was best described by means of a numerical 
function (pointwise) determined from the measurement 
of the 208pb sample (see Ref. [H] for details). The un- 
certainties quoted for the energy of each resonance are 
only the statistical errors obtained from the fits of the 
capture data performed with SAMMY. 



TABLE II: Resonance parameters derived from the R-matrix analysis of 
the 2"®Pb(n,7) data. 



Eo 


I 


J 






Fn 


AF„ 






(eV) 






(meV) 


(%) 


(meV) 


(%) 


(meV) 


(%) 


3357.93(0.04) 


1 


3/2 


78.1 


3 


235 




117 


2 


10865.0(0.4) 


1 


3/2 


64.9 


9 


44.1 


8 


52.5 


6 


11296.0(0.5) 


(1) 


(1/2) 


455 




44.6 


7 


40.6 


7 


14220.0(0.6) 


1 


(1/2) 


152 


6 


1560 




139 


5 


16428.0(0.4) 





1/2 


2268 


9 


936 


5 


662 


5 


19744.0(1.3) 


1 


(1/2) 


156 


7 


2581 




147 


6 


19809.0(0.9) 


1 


(3/2) 


295 




71.6 


8 


115 


6 


21885.0(0.9) 


1 


3/2 


121 


6 


875 




212 


5 


25112.0(0.9) 


1 


3/2 


438 


9 


326 


8 


374 


6 


25428(5) 


1 


1/2 


254 


7 


48901 




253 


7 


36200(6) 


1 


1/2 


312 


14 


35700 




309 


13 


37480.0(1.9) 


1 


(3/2) 


151 


15 


890 




258 


13 


39028(2) 


1 


(1/2) 


346 




93.0 


36 


73.3 


28 


40647(2) 


1 


(1/2) 


163 


23 


884 




138 


19 


42083.0(1.7) 


1 


(3/2) 


419 


21 


1419 


91 


647 


26 


47534(2) 


(1) 


(1/2) 


184 


34 


1000 




155 


29 


59233.0(0.2) 


(2) 


(3/2) 


322 


16 


1000 




487 


12 


63976(3) 


(2) 


5/2 


151 


17 


1110 




400 


15 


65990(10) 





1/2 


1186 


9 


82200 




1169 


9 


66590(6) 


1 


3/2 


198 


19 


9530 




387 


19 


70352(7) 


1 


1/2 


163 


34 


10780 




161 


34 


80388(4) 


2 


3/2 


1490 


8 


7005 




2457 


6 


83699(6) 


(2) 


(3/2) 


351 


16 


8000 




673 


15 


88509(6) 


2 


5/2 


375 


13 


7996 




1076 


12 


91740(4) 


(1) 


(3/2) 


298 


25 


1000 




460 


19 


92620(13) 





1/2 


991 


15 


32000 




961 


15 


93561(6) 


2 


3/2 


125 


37 


7001 




246 


37 


94743(7) 


2 


(3/2) 


241 


20 


7000 




465 


20 


101220(7) 


2 


(5/2) 


119 


26 


8000 




351 


25 


114380(5) 


1 


(3/2) 


655 


24 


2500 




1037 


19 


114602(6) 


2 


(5/2) 


366 


19 


5600 




1030 


18 


118100(6) 


2 


(5/2) 


390 


16 


5100 




1087 


15 


124753(47) 


1 


3/2 


2972 


9 


300000 




5886 


9 


125312(7) 


2 


(3/2) 


2783 


10 


21005 




4915 


9 


126138(38) 


(1) 


(3/2) 


319 


32 


100000 




635 


32 


140570(23) 


2 


3/2 


1387 


11 


103000 




2736 


11 


145201(6) 


(2) 


(3/2) 


518 


30 


3100 




888 


26 


146419(24) 





1/2 


6092 


8 


176000 




5888 


8 



7 



150880(7) 


(1) 


(1/2) 


554 


48 


4400 


492 


43 


151290(13) 


2 


5/2 


457 


23 


19000 


1340 


22 


191217(48) 


(1) 


(1/2) 


767 


28 


96977 


761 


27 


196990(37) 


1 


1/2 


584 


45 


64000 


579 


44 


198618(34) 


2 


3/2 


2730 


10 


132108 


5350 


10 


274630(22) 


1 


(1/2) 


514 


65 


32000 


506 


64 


276984(49) 


2 


3/2 


2481 


13 


112000 


4854 


13 


313400(18) 


2 


(3/2) 


1020 


32 


22000 


1950 


31 


314340(84) 


2 


5/2 


964 


24 


179000 


2875 


24 


356098(22) 


2 


(5/2) 


676 


35 


31000 


1985 


35 


357465(87) 


2 


3/2 


1998 


24 


455000 


3979 


24 


406200(55) 


2 


5/2 


656 


51 


102000 


1955 


51 


407200(41) 


2 


3/2 


2906 


24 


71000 


5583 


23 


416370(127) 


2 


5/2 


2722 


16 


307000 


8096 


16 


433340(32) 


2 


(5/2) 


4122 


17 


47000 


11368 


16 


434604(37) 


(2) 


(3/2) 


4695 


23 


58000 


8687 


21 


443412(13) 


(2) 


(5/2) 


2375 


21 


14000 


6092 


18 


466320(49) 


(1) 


(3/2) 


5413 


15 


90000 


10211 


14 


469080(76) 


2 


3/2 


3222 


19 


161000 


6317 


19 


471789(28) 


(3) 


(5/2) 


792 


36 


41000 


2330 


35 


476310(172) 





1/2 


5252 


18 


374000 


5180 


17 


510690(51) 


(2) 


(3/2) 


3123 


18 


86000 


6026 


18 


572245(181) 


2 


5/2 


3838 


13 


793194 


11460 


13 



" Capture kernel Kr = 5r^r„/r, with g = J + 1/2. 



A. Comparison to previous work 

The radiative neutron capture cross section of ^°^Pb 
has been measured at ORNL [i3,[ll[23|, at RPI [H and 
at IRMM 26] . As representative examples of these mea- 
surements we consider in this section two measurements 
made at ORELA [3, [3 ' ^ more complete analysis [l^] 
of the ORELA capture data made in combination 
with transmission data [28l | and the recent experiment 
made at IRMM [2^1 ■ In order to compare these four data 
sets with the present results (Table HI)), the ratio of the 
capture kernels are shown in Fig. [3l 

The values reported in Ref. ^J^] show a relatively good 
agreement with our results, except for the first two res- 
onances at 3.3 keV and 14.25 keV, which are lower by 
~50% (see Fig. [3]). However, these two resonances and 
the resonance at 16.428 keV are important because of 
their dominant contribution to the MACS in the energy 
range between 5 keV and 20 keV. It is difficult to deter- 
mine the source of discrepancy, thus no correlation has 
been found between the discrepancies and the spins of the 
resonances. The latter could probably help to determine 
if there is any effect related to the angular distribution 
of the prompt capture 7-rays or to the WF used in the 
previous measurement. 

In the second measurement at ORELA 141] the discrep- 
ancies versus our present results are smaller (see Fig. [3|) , 
but the capture kernels are systematically larger, on av- 
erage 20±5% higher. This could probably reflect that the 
WF used in Ref. [3] is overweighting the relatively hard 
pulse height spectrum of ^"^Pb. Indeed, similar discrep- 
ancies have been found in the past for ^^Fe [2^, where the 
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FIG. 3: (Color online) Ratio between the capture kernels re- 
ported in Ref. [l|] (top), Ref. [Ti] (2nd), Ref. ;27i (3rd) and 
Ref. [2^ (bottom) and the kernels determined here. 



pulse height spectrum is also considerably harder than 
that of the ^^''Au sample used for yield normalization. 

The posterior analysis [l^l of the ORELA capture 
data [3] in combination with transmission [2^ shows, 
on average, better agreement with the capture areas re- 
ported here (see Fig. [3]). 

Finally, the results reported in the measurement at 
IRMM [2y| show the best agreement with the capture 
kernels of n_TOF (see Fig. At E„ < 40 keV both 
measurements agree within a few percent. At higher en- 
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ergy the fluctuations are larger, but the agreement is still 
good within the quoted error bars. 

As an illustrative example, the capture yield measured 
at n_TOF for the first resonance at 3.3 keV is compared in 
the top panel of Fig. Inversus the yield calculated from the 
resonance parameters reported in Refs. [13, IH, [2^ . Ob- 
viously, the IRMM and n_TOF results show good agree- 
ment in both the capture area and the resonance energy. 



2 

(U 

>0.15 



0.1 



0.05 



T77/S Work 




\RMM, 2007 




- ORELA, 1979 / \ 




Mughabghab 2006 /7\ In 






*'.*i-'VWA'^,j 


3.35 3.36 


3.37 



En (keV) 



(U 

0.03 



0.02 



0.01 




15 



FIG. 4: (Color online) (Top figure) The bold red line repre- 
sents an R-matrix fit to our experimental capture yield start- 
ing from the initial parameters (solid green line) in Ref. [23 |. 
The dashed and dott-dashed curves correspond to the cap- 
ture yields determined in Ref. [2^ and Ref. [l^ . respectively. 
(Bottom figure) The fitted capture yield in the 10-30 keV en- 
ergy range (thin red line). 



B. Maxwellian averaged capture cross section 

The Maxwellian averaged cross section (MACS) was 
determined using the SAMMY code in the range of ther- 
mal energies relevant for stellar nucleosynthesis, i.e. from 
kT=5 keV up to fcT=50 keV. As discussed in the previ- 
ous section, our results agree best with the values re- 



ported in Ref. [26[. The latter data set seems also to 
be the most complete in terms of number of analyzed 
resonances, with about 283 levels. Therefore our results 
were complemented with resonances from Ref. ^265 in or- 
der to avoid any discrepancy due to resonances missing 
in Table Hn The contribution of these supplementary res- 
onances to the MACS is < 0.1% at fcT = 5 keV and 6% 
at kT = 25 keV. The fact that this correction starts to 
be significant towards kT > 25 keV is not relevant for 
the study of the nucleosynthesis of ^°^Pb. Indeed, as it is 
discussed below in Sec.|Vl ^"^Pb is mostly synthesized be- 
tween the He-shell flashes of the asymptotic giant branch 
stars. These intervals between pulses provide about 95% 
of the neutron exposure via the ^^C(a,n)^^0 reaction, 
which operates at a thermal energy of kT = 8 keV. At 
this stellar temperature less than 0.5% of the MACS is 
due to the supplemented resonances. 



This work 
IRMM 2007 
Mughabghab'06 
Bao et al. 
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FIG. 5: (Color online) Maxwellian averaged (n, 7) cross 
sections for ^'^'^Pb from the resonance parameters of this 
work (bold red) compared to the IRMM measurement [26l ] 
(dashed) , to the recommended data of Ref. [s^ (grey) and to 
the compiled data of Ref. (solid green). 

The uncertainties shown in Fig. [5] are only statistical. 
The systematic uncertainties of the MACS quoted in Ta- 
ble |TTT] include all contributions discussed in Sec. IIII Dl 



TABLE III: Maxwellian averaged cross section for ^""^Pb. 



Thermal energy kT 


MACS 


C" Stat 


^ sys 


(keV) 


(mbarn) 


(%) 


(%) 


5 


21.3 


1.8 


5 


8 


20.4 


1.8 


3 


10 


19.4 


1.9 


3 


12 


18.4 


2.0 


3 


15 


17.1 


2.1 


3 


20 


15.6 


2.2 


3 


25 


14.7 


2.3 


3 


30 


14.2 


2.3 


4 


40 


13.5 


2.2 


4 


50 


12.8 


2.1 


4 



Assuming systematic uncertainties of 4% and 10% for 
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3/2~ and 5/2~ resonances, respectively, the final uncer- 
tainties are completely dominated by the 4% uncertainty 
of the 3/2~ resonances. A change of 10% in the cross 
section of the fewer 5/2+ resonances has a negligible in- 
fluence on the MACS at kT = 5 keV, it contributes only 
0.5% at kT — 25 keV and increases linearly up to 1% at 
kT = 50 keV. An effect of 10% in the capture yield of 
the 3/2+ resonances makes only a 1% difference in the 
MACS at kT = 25 keV and it becomes also negligible 
towards lower stellar temperatures. The 3% systematic 
uncertainty of the experimental method itself originates 
from the PHWT, the neutron flux shape, and the use of 
the saturated resonance technique. 

In summary, the MACS of ^°^Pb can now be given 
with total uncertainties of 5% and 4% at the stellar tem- 
peratures corresponding to 5 keV and 25 keV thermal 
energies, respectively. This improvement with respect to 
the previously recommended values of Ref. [l^l becomes 
particularly important for determining the s-process con- 
tribution to the production of lead and bismuth in the 
Galaxy. 



THE s-PROCESS ABUNDANCE OF ^'^'^Pb AS 
A CONSTRAINT FOR THE U/TH CLOCK 



ation between observed and calculated abundances for 
s-process only isotopes This uncertainty is justified 
for ^°^Pb because its nucleosynthesis is dominated by 
the main component and it is only marg inally affected 
(~2%) by the strong component ^l|, l32l [33l [3^ . Fur- 
thermore, because of the much lower cross sections of 
•^"^Pb and ^°^Bi, the synthesis of ^°^Pb remains practi- 
cally unaffected by the a-recy cling after ^°^Bi 7]. This 
lends further confindence that the production of ^°^Pb, 
and hence its uncertainty, follows the same trend as the 
main s-process component. 
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The s-process production of 206p|-, ^^kes place in low 
mass asymptotic giant branch (agb) stars of low metallic- 
ity [sH, where about 95% of the neutron exposure is pro- 
vided by the ^^C(q!, n)-'^^0 reaction at a thermal energy 
of kT w 8 keV. At this stellar temperature the present 
MACS is about 20% lower and two times more accu- 
rate (see Fig. [5|) than the values from Ref. [13], which 
have been commonly used so far for stellar nucleosyn- 
thesis calculations. The additional neutron irradiation 
provided by the ^^Ne(a, n)^^Mg reaction at the higher 
thermal energy of kT ~ 23 keV during the He shell flash 
is rather weak. 

With the new MACS the s-process abundance of ^"^^Pb 
has been re-determined more accurately. A model calcu- 
lation was carried out for thermally pulsing AGB stars of 
1.5 and 3 Mq and a metallicity of [Fe/H] = -0.3. The 
abundance of ^"^Pb is well described by the average of 
the two stellar models, which represent the so-called main 
component [s^]. Since the contribution of ^°^Pb by the 
strong component is only 2%, the main component can 
be used to approximate the effective production of ^'^^Pb 
during Galactic chemical evolution (GCE) [3l|, [H, Is^ . 
This approach yields an s-process abundance of ^°^Pb, 
which represents 70(6)% of the solar abundance value 
iV206 = 0.601(47)/106Si [HI. The same calculation made 
with the older MACS recommended by Bao et al. [s^ 
yields 64%. The uncertainty on the calculated s-process 
abundance is mostly due to the uncertainty on the solar 
abundance of lead (7.8%) [1^. The contribution from the 
uncertainty on the MACS at 8 keV is less than 2%. Fi- 
nally, the contribution from the s-process model is ±3%. 
The latter corresponds to the mean root square devi- 



FIG. 6: Estimate of the radiogenic component of Pb using 
the Fowler's model with different nucleosynthetic assumptions 
(see labels in curves) and the r-process age = tu — 4.6 Gyr 
(vertical dashed line) derived from the age of the Universe 

tu [131. 

In order to estimate a constraint for the r-process 
abundance of ^°^Pb one needs to take into account its ra- 
diogenic contribution, Nf^^, due to the decay of ^^^U. As 
it is shown in the following, this component is relatively 
small but can not be neglected. Based on the schematic 
model of Fowler, which assumes an exponential decrease 
of the r-process yield during GCE 4] (supernova rate A = 
(0.43t)~^ Gyr~^) and using the current best estimates for 
the age of the Universe {tu = 13.7 ± 0.2 Gyr) [37*1, one 
obtains A^^oe ^ o.027(2)/10^Si (see Fig.lHand Table[IVl). 
This number, combined with our resuh for A^^oe^ yields 
an r-process residual, 

TV^oe ^ ATgoe _ _ ^ 0.153 ± 0.063. (6) 

The uncertainty in this result includes contributions of 
8.4% from Nf^^ (corresponding to the uncertainty on the 
solar abundance of ^ssy [35|)^ 7.8% from the total solar 
abundance of ^o^pb, A^|06^,|3^, and 8.6% from the de- 
termination of A^"^ as discussed above. This means that, 
apart from the uncertainties related with the simplified 
assumptions in the GCE model of Fowler, the r-process 
abundance can be reliably constrained between 16% and 
36% of the solar ^ospb. 

The r-process residuals derived here are consistent 
with r-process model calculations available in the liter- 
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TABLE IV: Radiogenic abundance of ^^^Pb, iV^^'"' (Si^lO''), 
derived from the model of Fowler and the age of the Universe 
(see Fig. [6)l . r-Process residuals obtained via eq. |6] 



GCE 
(Fig.® 



lO^'Si 



lO^'Si 



N? 



43% SN rate 

Sudden 

Uniform 



0.027(2) 
0.058(5) 
0.0161(14) 



0.15(6) 
0.12(6) 
0.16(6) 



26(10) 
20(10) 
27(10) 



TABLE V: Radiogenic abundance of ^''''Pb, N^°'' (Si=10''), 
derived from the model of Fowler and the age of the Universe. 
r-Process residuals obtained via eq. [S] 



GCE 



lO'-Si 



lO'^Si 



■207 /^r207/(W- 



43% SN rate 
90% SN rate 
Uniform 



0.150(13) 
0.08(7) 
0.047(4) 



0.003(73) 
0.073(72) 
0.106(72) 



0(11) 
11(11) 
16(11) 



ature, i.e., iV^°^ = 26.6% More recent calculations 
yield TV^oe values between 27% and 35% d. One can 
also derive hard limits for the r-process abundance, con- 
sidering the two extreme cases of sudden nucleosynthesis 
(A — ^ oo) and uniform nucleosynthesis (A — > 0). This 
yields constraints between 10% and 37% of solar ^"^Pb 
(see Table HyI). 

The situation is rather different for the correspond- 
ing ^°^Pb/^^^U ratio, which has been investigated as 
a potential clock in the past [38]. In this case, the 
s-process abundance of ^°^Pb was recently determined 
to be Nf^ = 77(8)% A similar calculation to 

that shown in Fig. M gives N^°'^ = 0.150(13) (see Ta- 
ble |V|. The latter value reflects the large relative ra- 
diogenic abundance of ^o^Pb, N^07/j^207 ^ 22%, due to 
the much shorter half-life of ^'^^U. From the total solar 
abundance of ^o^Pb ^ and the N^'^ and iV^o^ values 
quoted above, the r-process residual becomes N^'^'^ = 
0.003 ± 0.073, which means that iV^°^ can not be larger 
than 11% of the ^°^Pb abundance in the solar system, 
Ar207 ^ 0.665(52) [H (Table|V|. 

This result is in contrast with r-process model calcula- 
tions, which yield values between 22.7% and 25.3%, with 
a relative uncertainty of 15-20% [l|, The s-process 
abundances of 206,207pb ^re rather reliable and not very 
sensitive to details of the stellar models 0, ■ There- 
fore, this discrepancy indicates that r-process abun- 
dances might have been overestimated, possibly because 
the odd-even effect is not properly reproduced by the 
ETFSI-Q mass model implemented in the r-process cal- 



culations [U, Q. Indeed, one needs to increase the super- 
nova rate in the standard Fowler model from 43% up to 
90% (A = (0.90t)~^ Gyr^i) in order to achieve agree- 
ment between these r-process constraints and the latter 
r-process calculations [l|, Q . Obviously the less realistic 
uniform scenario would also provide agreement with the 
abundances from these r-process models (see Table |Vl . 

However the situation has been improved recently after 
more detailed r-process model calculations [i^l predicted 
a new A'^'^^ value, which is 35% lower than the previous 
one of Ref. Q. This yields Af207/^207 ^ ^^6.8%, which is 
substantially closer (considering an uncertainty of 20%) 
to the upper hmit of 11% derived here. In this case a good 
agreement would be found for a more reasonable increase 
of the supernova rate to 55% in the Fowler model. 

These constraints for the r-process abundances of 
206,207p|-j become relevant for the validation of r-process 
model calculations and hence, for the reliable interpreta- 
tion of actinide abundances observed in UMP stars and 
their use as cosmochronometers. 

The s-process aspects will be more rigorously investi- 
gated in a comprehensive study of the Pb/Bi region [4]| . 
where the role of stellar modeling and GCE will be dis- 
cussed with a complete set of new cross sections for the 
involved isotopes, includin g th e present data for ^°^Pb, 
and recent results for 204pb [2|, 207pb [l| and ^"SBi [ll|- 



VI. 



SUMMARY 



The neutron capture cross section of ^"^Pb as a func- 
tion of the neutron energy has been measured with high 
resolution at the CERN n_TOF installation using two 
CeDe detectors. Capture widths and/or radiative kernels 
could be determined for 131 resonances in the neutron 
energy interval from 3 keV up to 620 keV. Systematic 
uncertainties of 3%, 5%, and <10% were obtained for 
resonances with spin-parities of 1/2='=, 3/2~ and 5/2+, 
respectively. The Maxwellian averaged cross sections 
were found to be significantly smaller by 10% to 20% 
compared to values reported earlier [30|j resulting in a 
correspondingly enhanced s-process production of ^"^^Pb. 
First calculations with a standard AGB model yield an s- 
process component of 70(6)% for the ^°^Pb abundance. 
Combined with an estimate of the radiogenic produc- 
tion of ^"^Pb, the r-process abundance is constrained be- 
tween 16% and 36% of the solar ^°^Pb abundance, well 
in agreement with r-process model calculations reported 
in the literature ^, 's'l . A similar analysis for ^°^Pb shows 
agreement only with most recent r-process model calcu- 
lations [H. 
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